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In order to later find explicit analytic solutions, we investigate the singularity 
structure of a fundamental model of nonlinear optics, the four-wave mixing 
model in one space variable z. This structure is quite similar, and this is not 
a surprise, to that of the cubic complex Ginzburg-Landau equation. The main 
\J^} ■ result is that, in order to be single valued, time-dependent solutions should de- 

■ pend on the space-time coordinates through the reduced variable £ = ^fze~ t l r , 

qq ' in which r is the relaxation time. 

O 
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1. The four- wave mixing model 

Dynamic holography relies on the ability of the refractive index of a non- 
linear media to locally change under the action of light, i.e. under the so- 
called photorefractive response. The most conventional scheme for dynamic 
holographic applications consists in the formation of a fringe interference 
pattern between the interacting waves (the light grating) , which creates the 
refractive index grating. During the wave interaction, these are the same 
waves which create the grating pattern and diffract on this grating. There 
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exist many applications based on dynamic holography with photorefractive 
crystals 1 . 

In the four-wave mixing model, the grating, i.e. the modulation of the 
refractive index induced by the light interference pattern, is created by two 
pairs of co-propagating waves, 1-2 and 3-4. Let us denote A\ , A 2 , A 3 , A 4 
the complex amplitudes of the four waves, Ae the complex amplitude of 
the index grating, r the characteristic time of the evolution, I the total 
intensity of light (a real positive constant), 7 the photorefractive coupling 
(a complex constant). The four- wave mixing model is the system of five 
complex partial differential equations obeyed by these five complex ampli- 
tudes 2 (bar denotes complex conjugation) 

d z A x = -iAeA 2 , d z A 2 = iAeAi, d z A 3 = -iAeA 4 , d z A 4 = iAeA 3l 

d t + -)Ae = ^- (A{A 2 + A 3 A 4 ) . (1) 

The correspondence of notation with Ref. 3 is 

{A 1 ,A 2 ,A 3 ,A 4 ,Q, 1 ) of Ref. = (A u A 4 ,A 3 ,A 2 ,-iAs,ij) hcic . (2) 

At present time, no analytic solution is known to this system, except a 
stationary solution in which the grating amplitude has a pulse profile 1 , 

Ae = Kscchk(z — z ), (K, k) — real constants. (3) 

The purpose of this short article is to investigate the structure of the sin- 
gularities of the time-dependent solutions (dt 7^ 0) in order to find, in a 
forthcoming paper, closed form analytic solutions by making suitable as- 
sumptions dictated by the singularity structure. 

The first four equations (1), which do not depend on dt, admit six 
quadratic first integrals. Two of them, 

|Ai| 2 + \A 2 \ 2 = fl{t), |A 3 | 2 + |A 4 | 2 - /f (i), (4) 

allow us to represent the four complex amplitudes Aj(z,t) with two real 
functions /1, f 2 of t and six real functions 9\, 9 2l tpi, ip 2 , ip 3l ip 4 of (z, t), 

A x = fi sin 9 ie 1 ^ ,A 2 = h cos B x e 1 ^ , A 3 = f 2 cos 6 2 e 1 ^ , A 4 = f 2 sin 9 2 e l ^{5) 

The four other quadratic first integrals are pairwise complex conjugate 




A 2 A 4 = f 3 (t), AxAs + A 2 A 4 = h(t), 
A 2 A 3 = f 4 {t), A X A 4 - A 2 A 3 = f 4 (t), { ' 

and constrained by one relation, 

-fill + l/ 3 | 2 + |/ 4 | 2 = 0. (7) 
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2. Local singularity analysis of the FWM model 

The dependent variables Aj, Ae generically present movable singularities, 
i.e. singularities whose location in the complex domain depends on the ini- 
tial conditions. The study of the behaviour of Aj,Ae near these movable 
singularities is a prerequisite to the possible obtention of closed form ana- 
lytic solutions to the nonlinear system (1), so let us first perform it. Such a 
study is made by performing the successive steps of the so-called Painleve 
test, a procedure explained in detail in Rcf. 4 ' 5 . 

The movable singularities of a PDE lay on a manifold represented by 
the equation 6 



ip(z,t) -ip o = 0, 



(8) 



in which ip is an arbitrary function of the independent variables, and ipo 
an arbitrary movable constant. In order to express the behaviour of the 
ten dependent variables A, , A, , Ae, Ae, it is convenient to introduce an 
expansion variable x(z, t) 4 which vanishes as ip — ipo an d which is defined 
by its gradient 



i S 2 

Xz = 1 + , Xt 
with the cross-derivative condition 



-c + c zX --(cs + c zz )x 2 . 



X = Sf ~t~ Czzz *2-CzS ~\~ CS Z 



0. 



(9) 



(10) 



Since they will be needed later, we give the dependence of x, S, C on tp, 

x ■ 



Vz 



2<p z 



Vz 



<Pz 



<Pzz 

2 \ V 



2 



c 



,^z^Q, (11) 

<Pt/<Pz- (12) 



2.1. Leading behaviour 

The first step is to determine all possible families of noncharacteristic 
(i.e. (p x (ft ^ 0) movable singularities, i.e. all the leading behaviours 



X 



: 



A k ~ a kX Pk , A k ~ b kX Pk , a k b k ^ 0, k = 1, 2, 3, 4, 
Ae ~ q X a , Ae ~ ■, qor + 0, 



(13) 



in which the exponents p k , P k , a, (3 are not all positive integers. If one as- 
sumes that the two terms A1A4 and A3A2 have the same singularity order, 
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the ten exponents pk, Pk, ct, (5 obey the twelve linear equations, 



Pi - 1 = p 2 + a, P 2 - 1 = Pi + a, 
Pa - 1 = p 3 + a, P 3 - 1 = Pa + a, 
P 1 -l = P 2 +0, P2-l=Pl+0, 
P 4 -l = P 3 + (3, p 3 -l=P4 + /3, 
a - 1 = pi + P 2 = Pa + P3, 
J3-l = P 1 +p 2 =P 4 +p 3 , 

and the ten coefficients a,k,bk,qo, ro obey the ten nonlinear equations 



(14) 



piai = -q a 2 , p 2 a 2 = -r ai, p 4 a 4 = -qoa 3 , p 3 a 3 = -r a 4 , 
Pi 61 = r b 2 , P 2 b 2 = q bi, P 4 b 4 = r b 3 , P 3 b 3 = q b 4 , 
Caq = -^y/I ){a\b 2 + a 4 b 3 ), 
k C(3r = -(7//o)(6ia2 + b 4 a 3 ). 



(15) 



Therefore the five squared moduli \Ak\ 2 , |Ae| 2 behave like double poles, 



Pk + Pk = -2, k = l, 2,3,4, a + = -2, 

and the linear system (14) is solved as 

f Pi = -1 + si + <5, Pi = -1 - si - <5, 
p 2 = -l + si - 5, P 2 = -1 - si + (5, 
p 4 = — 1 + s 2 + 5, P4 = — 1 — s 2 — 5, 
p 3 = -1 + s 2 -S, P 3 = -1 - s 2 + 5, 
la = -1 + 2(5, /?= -1-25, 

in which <5, s\, s 2 are to be determined by the nonlinear system (15). 
The first set of relations 

-Qor Q = P1P2 = P4P 3 = P1P2 = P4P3, 

imply 

si =0, s 2 = 0, -g »*o = 1 - <5 2 . 

The second set of relations 

01 a 4 b 2 b 3 iq 1 + 5 



(16) 



(17) 



&2 



03 



61 b 4 1 - <5 



zr 



is solved as 



ai = ai 2 \, b 2 = -612A, 04 = a 34 A, 63 = 634A, iqo = (1 - <5)A 2 , 
a 2 = ai 2 A _1 , 61 = &i 2 A _1 , a 3 = -a 34 X^ x , b 4 = b 34 \~ x , ir = (1 + <5)A" 



(18) 
(19) 
(20) 

^1) 



in which a\ 2 , a 34 , b\ 2 , b 34 , A are nonzero functions of (z, t) to be determined. 
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Finally, the two remaining equations 



A12&12 " 034&34 + — (1 - S){\ - 25) = 0, 
ai2&i2 - a 34 fo 3 4 - + 5)(1 + 25) = 0, 



first define 5 as a root of the second degree equation 
7(1 + 5)(1 + 25) + 7(1 - 5){l - 25) 



0, 



(22) 



(23) 



then put one constraint among a\2, 612, 034, 634- 

The equation (23) which defines 5 only depends on the argument of 7, 

(2 cos g)5 2 + (3i sin g)5 + cosg = 0, ,g = arg7. (24) 

When the photorefractive coupling constant 7 is purely imaginary, the ex- 
ponent 5 vanishes, otherwise it can take two purely imaginary values 



. -3 sing ± ^9 sin 2 g + 8 cos 2 g 

5 = 1 - . (25) 

4 cos g 

To conclude, there generically exist two families (13) of movable singu- 
larities, defined by the equations (17), (19), (21), (22), i.e. 

5 = one the two roots of Eq. (24), 
Pi =pa = P 2 = P3 = -1 + 5, 
Pi = Pa = Pi = P3 = -1 - 5, 
(Zi = iVA P12 cosh ji, 62 = — iVA cosh /x, 
a4 = iVA p 3 4 sinh^i, 63= iVA p^smhfi, 
a 2 = N\~ 1 pi2 cosh/x, bi — N\~ 1 p^ 2 1 cosh [i, (26) 
a 3 = — 7VA _1 |J34 sinhp, 64= iVA -1 ^ 1 sinh^i, 
i<7o = (1-5)A 2 , zr = (l + r5)A- 2 , 

^-^(l-^Xl-tf 8 ), 
G = (1 + <f)(l + 25) 7 = -(1 - - 25)7- 
and they depend on four arbitrary complex functions A, fi,pi2,P34 of (z, t). 



2.2. Fuchs indices 

The leading behaviour (13) is the first term of a Laurent series 



X^0: { 



A k = X Pk J2 a k,iX j , &£ = X a J2 9jX j , 



+OO 



j=o j=0 



(27) 
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and the indices j at which arbitrary coefficients enter this expansion are 
computed as follows. If one symbolically denotes E{A k , A k , Ae, Ae) = 
the ten equations (1), one builds the linearized system of (1), 

Um E{A k + (C k ,A k + (C k , Ae + (D, Ae + (D) - E(A k ,A k , Ae, Ae) = Q 
i.e. 

d z Ci + iAeC 2 + iA 2 D = 0, ... 

(dt + D - ^ (3 2 Ci + A 1 C 2 + A 3 C 4 + A 4 C 3 ) = 0. (28) 

At the point (13), this ten-dimensional linear system displays near \ — a 
singularity which has the Fuchsian type 7 , therefore it admits a solution 

n / C k ~ c kX Pk+j , C k ~ d kX Pk+j , k = 1,2,3,4, 

The condition that this solution be nonidcntically zero results in the van- 
ishing of a tenth order determinant whose roots j, called Fuchs indices, take 
the ten values 

j = -1,0,0,0,0,2,2,2, i±^jgg (30) 

In the two-wave mixing case, the four roots 0,0,2,2 are to be removed 
from this list. The necessary condition, required by the Painleve test, that 
all indices be integer is violated when 6 is nonzero, therefore the FWM 
model fails the test in this case. This failure is quite similar to what occurs 
in the complex cubic Ginzburg-Landau equation 8 , and conclusions similar 
to those of Ref. 9 can probably be drawn. 

Let us remark that the Fuchs index has multiplicity four, in agreement 
with the number of arbitrary functions involved in the leading behaviour, 
see Eq. (26). 



2.3. Conditions at the Fuchs indices 

The coefficients a k j,b k j,qj,Tj of the Laurent series (27) are computed for 
j > 1 by solving a linear system. Whenever j reaches an integer Fuchs 
index, an obstruction may occur, resulting in the introduction of a movable 
logarithmic branching, and some necessary conditions need to be satisfied 
in order to avoid it. In the generic case 5^0, these obstructions may only 
arise at j = 2, value of a triple Fuchs index. In the nongeneric case 8 = 0, 
in addition to the quadruple index j = 2, one must also check the simple 
index j = 3. The results are as follows. 
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At j = 2, whatever be 6, a movable logarithm exists, unless the following 
necessary condition is satisfied, 

Q 2 = a{C t +CC z -2aC) = Q 1 o = -. (31) 

r 

Since a = is excluded, the second factor defines a precise dependence of 
C on (z,t), hence a similar dependence for S,x,f arid ultimately for the 
five complex amplitudes via their Laurent expansion (27). Let us find this 
dependence explicitly. 

The general solution of (31) is provided by the method of characteristics 
and is defined implicitly by the relation 

2az = C + F(e~ 2at C), (32) 

in which F is an arbitrary function of one variable. 

The invariant S is then defined by the first order linear PDE (10), which 
reads 

S(x, t) = E(C, t), D = l + e- 2at F', 

(v \ F" 2 p m 

S t + 2aC£ c + 4a-j+ 24a 3 e ~ 6at — - 8a 3 e~ 4at ^ - 0,(33) 

and admits the general solution 

S = c e- iat D- 2 - Aa 2 e^ at D^F"' + %a 2 e~ &at D' A F"\ (34) 

in which Co is an arbitrary constant. 

We have not been able to prove whether the function F is a gauge 
which can be arbitrarily chosen (e.g. F = 0) or whether it is essential. If 
it is essential, there could exist much more intricate solutions than those 
which we now outline. 

When F is arbitrary, we have not yet succeeded to integrate the system 
(9) for x(z, t), hence to compute the induced dependence of Aj, Ae on (z, t). 
However, if one restricts to a constant value for F, which can then be chosen 
equal to zero, the equation (32) for ip integrates as 

F = 0, ip = $(£), £ = ^ze- a \ (35) 

in which $ is an arbitrary function, and it is straightforward to check that 
the Laurent expansion (27) defines the reduction (z,t) — > £ = y/ze~ at . As 
opposed to the stationary reduction dt = 0, this reduction is noncharactcr- 
istic (i.e. it does not lower the differential order) it is defined with a reduced 
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variable of the factorized type £ = f(z)g(t), 

' i = z 1 ' 2 e- a \ o = 1/t, 
Ae(z,t) = e- o *- <w *«- 1 / 2 Ae r (0 ) Ae(«,t) = e -° t+ih '*«- 1 / 2 Ae r (0 ) 
Aj(z,t) = e-^-^A^), A^t) = e- at+ ^/ 2 A jtr (0, j = 1,4, 
A,-(s,t) - e- at+ -*/%, r (0, A,-(z,t) = e" *-*-*/ 2 ^^, J = 2,3, 
— Ai r = -2iAs r A 2r , —A 2r = 2iAe r Ai r , ( 36 ) 

d- - d 

-jzM,r = -2iAe r A^ r , ~rz A^. r = 2iAe r A 3 . r , 

d 

— Ae r + iuT^Asr = — — (A lr A 2r + A 3r A 4r ) , 

and it depends on the additional real parameter u>. 

To summarize the information obtained from the Fuchs index 2, pro- 
vided F can be arbitrarily chosen, the four-wave mixing model (1) (as well 
as the two- wave mixing model) admits no single valued dynamical solution 
other than the possible solutions of the reduction (z, £)—>£ = y/ze~ at . 

When S ^ 0, this ends the Painleve test. When 5 = 0, the Fuchs index 
3 is found to be free of movable logarithms. 

Finally, the information provided by the Painleve test is the following. 

(1) Whatever be 7, and under the mild restriction that F can be arbitrarily 
chosen, no single valued dynamical solution of (1) exists other than the 
possible solutions of the reduction (z, t) — > £ = y/ze~ at . 

(2) For this reduction, two cases must be distinguished. When 9?(7) 7^ 
0, the ten-dimensional ODE system (36) possesses at most an eight- 
parameter single valued solution. When 3?(7) = 0, the system (36) 
passes the Painleve test therefore it may admit a ten-parameter single 
valued solution, which then would be its general solution. Finding a 
Lax pair in this case would considerably help to perform the explicit 
integration. 

Similar conclusions apply to the six-dimensional two-wave mixing: the 
only possibility for a single valued solution is the reduction (36) (with A3 = 
A4 = 0), this solution depending on four movable constants for 3J(7) 7^ 0, 
and six movable constants for $£(7) = 0. 

3. Conclusion 

The present study proves, like for the complex Ginzburg-Landau equation, 
that physically relevant analytic solutions quite certainly exist for the four- 
wave mixing model. The present counting of the possible arbitrary constants 
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in the solutions displays, as expected, the crucial role of the photorefractive 
complex constant. Explicit solutions based on the present study will be 
presented elsewhere. 

The results can be used for predicting new phenomena in optical self- 
diffraction of waves in photorefractive media which use the non-local re- 
sponse as well as for optimization of optical dynamic holographic settings. 
Among these applications let us quote: (i) the formation of a localized 
grating to increase optical information density; (ii) the methods of all- 
optical control of output wave characteristics versus input beam intensi- 
ties and phases; (iii) the optimization of the parameters of optical phase- 
conjugation; (iv) the use of the new holographic topographic technique to 
material parameter characterization. Just like similar processes of nonlinear 
self-action of waves arise in models of optical networks, optical information 
processing, quantum information processing, interacted neural chains, then 
various other problems of nonlinear wave interaction can become the sub- 
ject of further independent research. 
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